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Abstract. We consider a simple stochastic model of a urban rental housing market, 

in which the interaction of tenants and landlords induces rent fluctuations. We 

, /^ i . simulate the model numerically and measure the equilibrium rent distribution, which 

^^ I is found to be close to a lognormal law. We also study the influence of the density 

pH ' of agents (or equivalently, the vacancy rate) on the rent distribution. A simplified 

version of the model, amenable to analytical treatment, is studied and leads to a 
lognormal distribution of rents. The predicted equilibrium value agrees quantitatively 
with numerical simulations, while a qualitative agreement is obtained for the standard 
deviation. The connection with non-equilibrium statistical physics models like ratchets 
is also emphasized. 



PACS numbers: 89.75.Fb, 02.50.Ey, 89.65.Lm 



^ ! 1. Introduction 



The emergence of a well-defined distribution of prices in a market involving a large 
number of interacting agents is a typical example of collective phenomenon in complex 



H ■ systems, for which concepts and tools borrowed from non-equilibrium statistical physics 

" " could prove useful. Theoretical studies of price fluctuations have been reported both 

in the economic [1] and in the physics literature, mostly in the framework of stock 
market models, where price fluctuations are large, with often fat-tailed distributions 
[21 [3l H]. Other markets, like the rental housing market or the labor market, seem to 
have attracted less attention from the physics community (see, however, [5l|6]), though 
their potential interest from a modeling point of view is certainly high. The economic 
literature on housing market is indeed in close relationship with the literature on labor 
market, which inspired some of the first models in the domain [3 Ej- These works study 
processes of search and matching between tenants (or persons wishing to purchase a 
house on the residential property market) and housing lots and their owners. The 
importance of the vacancy rate for the turnover between tenants and its influence on 
price formation have been underlined. Several works have studied price dispersion, and 
the conditions needed to model it on the housing market [9l |TOl [11], or corresponding 
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wage dispersion on the labor market [121 [T3] and general price dispersion for consumer 
goods [H]. The phenomenon of price dispersion is also documented in the empirical 
literature, be it for housing market [IHl HE] or prices of products in stores for instance 

The mechanisms at the origin of price dispersion in economic models are often 
different from that responsible for fluctuations in physical systems, namely thermal 
fluctuations. In most economic models, there is an exogenous heterogeneity in agents' 
behavior rules [11]: firms producing goods can have different costs of production, goods 
can have different characteristics [10], or consumers may have different search costs or 
behaviors [H]. The literature on search and matching theory focuses on the formation 
of relationships [HI [H [121 ^S\, between buyers and sellers for instance. There price 
dispersion usually comes from a difference in the information agents have on the market, 
or a difference in access to opportunities (e.g., a potential tenant is presented only a few 
fiat offers among many). 

Departing from this standard economic framework, we propose a parcimonious 
model of urban rental housing market, where the behavior of landlords and tenants is 
probabilistic, and driven by very limited information. Indeed, as explained in detail 
below, agents have in this model a very limited knowledge of the market, at odds with 
the "optimizing agent" framework often used in the search and matching literature. For 
instance, tenants only know their current rent, and the rent in a randomly chosen flat 
where they consider moving. Another important difference with the "optimizing agent" 
framework is our use of probabilistic rules for moves. We study both numerically and 
analytically this simple stochastic model, in order to characterize the rent distribution, 
and to determine the role of the density of agents (considered to be an external control 
parameter) on the rent formation mechanism. From a statistical physics viewpoint, our 
model also presents some similarities with non-equilibrium systems like ratchet models 
driven by dichotomous Markov noise [I9], suggesting interesting connections between 
the two fields. 

2. A simple model of rental housing market 

The model simulates the evolution of the distribution of rents in a city. Note that we 
focus on a rental housing market, and do not treat the question of the relation between 
rents and housing prices. In the model, the virtual city is composed of A'^; flats labeled 
by i = l,...,Ni, which are supposed to have identical characteristics, except for their 
rent, and A^ < A"; agents, with Na/Ni = p the density of agents. The flats can be in 
only two states: either occupied by an agent, or empty. When a flat is occupied, there 
is a probability for the tenant residing in it to move out, and when a flat is empty, there 
is a probability for a tenant to move in. 

In addition, the rent increases when a flat is occupied, as the rent is revised 
periodically by the owner. When a flat is empty, its rent decreases until some agent 
moves in. These ingredients give a simple dynamical framework to ensure a competition 
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on the housing market as in standard urban economics models [20], and are consistent 
with the fact that demand and supply drive the rent. From a physics point of view, 
this model is analogous to a lattice gas with long-range hopping, but is however more 
complicated due to the coupling with a dynamically evolving rent. 

The model is simulated in discrete time. At each time step, a flat is chosen 
at random in the city; this fiat is either occupied or empty. The other variable 
characterizing the fiat is its rent p, i.e., the rent that a tenant pays or would need to pay 
to occupy the fiat. Let us first consider the case where the chosen flat is occupied. In 
this case, there is a given probability vr^ that the rent is raised by the owner during the 
considered time step. If the rent is increased, it is multiplied by a given factor fr > 1. 
In a second substep, there is a probability tt^ that the tenant searches for another flat 
in the city with a lower rent. In this case, an empty flat, whose rent is denoted pe, is 
chosen at random in the city and the agent moves into this flat only if pe < p, the initial 
flat thus becoming empty. If Pe > p, the agent does not move, and the selected empty 
flat is rejected. The probability vr^ depends linearly on the rent p of the considered 
flat: TTj, = inin{p/ps,l), where Ps is a characteristic rent, which can be linked to the 
reservation rent in standard housing economics works. When the rent of the flat p is 
high, that is to say, close to or above Ps, the probability that the tenant looks for another 
flat is high. 

On the other hand, if the flat chosen at the beginning of the time step is empty, 
the rent is lowered by the owner with a probability tt/, in order to make the flat more 
attractive for potential tenants. At the end of the time step, the flat remains empty. 
This probability also depends on the rent p of the flat: vt/ = min(p/p/, 1), where pi is 
another characteristic rent, which could be related to rental provision costs in housing 
economics [10]. In case the rent is lowered, it is multiplied by a given factor /; < 1. Let 
us note that the behavior of landlords is also parcimonious. In order to optimize their 
profit, with a knowledge only of the state and rent of their flat, they raise the rent with 
a fixed probability when the flat is occupied (which can be interpreted as risk neutral 
behavior), but they are reluctant to decrease the rent when it is already too low. A 
summary of the dynamics of the model is given on Fig. [H 

3. Numerical simulations and results 

For practical reasons, and in order to facilitate comparison with the analytical model 
studied in Sec. HI it is convenient to discretize the decimal logarithm of the rent p. This 
discretization allows in particular for an easier determination of the rent distribution. 
We thus introduce a resolution a -C 1, such that log^gP = na, or equivalently p = lO'^", 
where n is a positive integer. For consistency, we also choose fr and fi in the form 
fr = 10""'° and // = 10~"'", with n^ and ni two given positive integers. 

The simulation is initialized as follows. Agents initially occupy the first Na flats 
i = 1, ...,Na, and the A'^; — Na remaining flats are empty. The rent of each flat is drawn 
at random from a normal distribution of given mean and variance. We find that the 
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Figure 1. Summary of a time step of the model. Note that tts = 1 ii p/ps > 1, and 
TT; = 1 if p/pi > 1 (see text). 
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Figure 2. Left panel: Evolution of the mean rent over time, with different initial rent 
distributions. Main figure: Gaussian initial distribution of the rent, with standard 
deviation 5 and mean 20, 40, . . ., 180. Inset: evolution of the mean rent for initial 
distributions with different shapes. Full line: Gaussian distribution of mean 90 and 
standard deviation 5; dashed line: uniform distribution between and 200; dotted 
line: Dirac peak at p=100. Parameters: Ni — 5000, Ps ^ Pi = 2000, tt^ = 0.03, 
fr = 10"" ~ 1.026, fi = 10-19" ~ 0.957 with a = 0.001, p = 0.7. Right panel: 
equilibrium distributions of the decimal logarithm a; of rent for the same initial 
distribution as on the inset of the left panel. All three distributions are identical 
to numerical accuracy. 



system reaches after some time an equilibrium rent that is independent of the initial 
rent distribution, as illustrated on the left panel of Fig. [2l 

To go beyond the average rent, and to characterize rent fluctuations, we measure 
the equilibrium rent distribution, which is obtained by averaging the instantaneous rent 
histogram over a long time lag, once equilibrium is reached. The resulting distribution 
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Figure 3. Left panel: averaged distribution <^{x) of the logarithm x of rent (for all 
flats, continuous line, for occupied flats only, dotted line, for vacant flats, dashed line) 
at equilibrium. Right panel: corresponding occupation rate r\{x) as a function of the 
logarithm x of rent (continuous line). Same parameters as Fig. [2l except for the right 
panel where ps = 1000 for the dotted line and p^ = 4000 for the dashed line. 



is presented in the right panel of Fig. [2l showing that the equihbrium distribution, and 
not only the mean rent, is the same independently of the initial distribution. Note that 
we actually measured the histogram of the logarithm of the rent, to better visualize the 
relative dispersion of rents. 

The left panel of Fig. [3] presents separately the distributions of the logarithm of 
rent for occupied and vacant flats, which shows that demand, corresponding to occupied 
flats, is centered around inferior rents than supply (vacant fiats), but both distributions 
intersect on a large range of rents, and build together the equilibrium rent distribution. 
This interplay of demand and supply leading to equilibrium rent provides a qualitative 
link between our work and standard economic theory on rent formation |21] . 

This feature of the model is also illustrated by another quantity of interest: the 
mean share of time during which a fiat is occupied, which characterizes the occupation 
of fiats at equilibrium (see right panel of Fig. [3l). This occupation rate illustrates the 
dynamics of the model: when rent is low, occupation rate is high, and vice versa. 
Occupation - corresponding to demand - is governing the rent of fiats in this model 
where rent is the only variable which distinguishes housing lots in the city. 

It is interesting to further explore the behavior of the model as a function of the 
density p of agents in the city, which is expected to have a significant impact on the 
equilibrium rent, since the density plays a key role in the demand. The evolution 
of the equilibrium rent distribution with p is shown on Fig. HI More precisely, the 
right panel of Fig. H] presents the rent distributions, while the left panel shows the 
distribution of the logarithm of the rent. The rent distributions are well fitted by 
lognormal laws (right panel of Fig. Hj), which simply appear as Gaussian laws when 
plotting the distribution of the logarithm of the rent. House rents have been often 
modelled as lognormal distributions in the housing economics literature following the 
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Figure 4. Equilibrium distributions (p{x) of the decimal logarithm of rent (left 
panel, continuous lines - dotted lines are Gaussian fits) and ip(j)) of rent (right panel, 
continuous lines - dotted lines are lognormal fits) for different values of the global 
density p. Same parameters as Fig. [21 except for the density: p = 0.1, 0.2, ..., 0.9 from 
left to right. 



founding work of [22], for example in mortgage default models [23]. Our model provides 
a parcimonious framework leading to this feature. 

Two main observations can be drawn from this study of the influence of density. 
The first one is the fact that the rent of housing increases when the density p is increased. 
This stylized fact is well-known in real cities [23], and is also reproduced by the standard 
urban economics model [20]. In the simple model presented here, this phenomenon can 
be interpreted as the fact that a higher density p yields a higher competition on the 
urban housing market, and thus higher rents. 

The second observation which can be made on Fig. H] is the decreasing width of 
the distribution of the logarithm of the rent, when the density p is increased. This 
can be linked to the multiplicative dynamics of the rent, according to the following 
interpretation. Let us suppose that rent p is such that p < Ps and p < pi. For high 
rents, a relative variation of the rent p (multiplication hj fr > 1 oi fi < 1) generates a 
higher absolute variation of the probability tt^ = p/ps that a tenant tries to move out of 
the flat (see the right panel of Fig. [3]) and of the probability ivi = p/pi that the rent of an 
empty flat decreases. These higher variations prevent the rent from increasing {ir,.) or 
decreasing (tt^) too much. Indeed, the reactions to rent variations are quicker for higher 
rents, and the rent is bound to stay within a smaller domain, in logarithmic scale. 



4. Analytical approach 



In order to better understand the numerical results obtained, we propose a simplified 
version of the model, that can be solved analytically. In this simplified model, we 
consider that the (decimal) logarithm x = log^^Q p of the rent p of housing describes a 
biased random walk on a one-dimensional lattice, with lattice spacing a, so that x = na. 
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Figure 5. Evolution of the mean x (left column) and of the standard deviation ax 
(right column) of the decimal logarithm x of rent as functions of the density p (top 
row) and the lattice spacing a (bottom row). Other parameters: same as in Fig. [5J 



This variable x evolves by integer jumps of rira = log^^Q fr in the positive direction, or 
nia = — log 5^0 fi i'^ ^^e negative direction, due to rent variations by the owner. In a 
mean-field spirit, the probability per unit time q for the rent to be raised should be 
proportional to the density p and to tt^, so that a first estimate would be g = iirp. 
Similarly, the probability q{x) for the rent to be lowered is proportional to 1 — p and to 

mix) =p{x)/pi = 10'' /pi, (1) 

leading to q{x) = (1 — p)ni{x). Note that only the probability to decrease the rent 
depends on the rent, through its logarithm x. 

To take into account the fact that a fiat remains occupied or empty during several 
successive time steps, so that the steps of the random walk should be correlated, we 
introduce the mean number Uo of rent increases during a period when a fiat is occupied, 
and the mean number Uy of rent decreases during a period when a fiat is vacant. We then 
consider that the variable x performs effective steps UoUra and nyUia (instead of Ura and 
riia), assumed to be statistically independent. In order to maintain the same average 
rate of rent increase and decrease, we divide the probabilities q and q{x) respectively by 
Ho and n^, yielding 



q = ^—, q{x) 



7ri{x){l- p) 



n. 



(2) 



Equating the mean increase rate qUoUr and decrease rate q{xeq)nyni of x, the equilibrium 
rent satisfies 



QiXeq) 



Q • 

muy 



(3) 
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Using Eqs. ([T]) and ([2]), we find the equilibrium rent Pcq = 10^'='': 

P=, ^±'^. (4) 

ni 1- p 

Note that the equihbrium rent is independent of the phenomenological parameters rio 
and n^. 

We now try to estimate the parameters Uo and riy. When a flat is occupied, at each 
step there is on average a probability approximately equal to ^vr^ = Pcq/i'^Ps) that the 
tenant moves to another flat, assuming that all rents are close to the equilibrium rent 
Peq- Thus the probability that a tenant has not changed fiats after n steps is 

-£)"• (^' 

and the probability that a tenant leaves the fiat after exactly n steps is /„ = Fn-i — Fn, 
yielding 

Hence, Uo can be estimated as the average number of steps a tenant stays in a given 
fiat, times the probability n^ that the rent is raised at each time step: Uo = i^r Yln'^^fn- 
Using Eqs. (j4]) and ([6]), this estimate leads to 

2psni{l-p) 
rio = . (7) 

PlUrp 

The computation of n^ follows the same scheme. The probability for a tenant to change 
fiats is Pcq/(2ps) at each time step, so that the probability for a vacant flat to welcome 
a tenant is 

^Vsii- p) 

(we assume in this simplified argument that the density is not too high, so that vr^ 
remains less than 1). The probability that the fiat remains vacant is 1 — 7r,„ at each 
time step, and (1 — vTu)" for n time steps. Following the same lines as for no-, we obtain 
n^ = 7ii{xeq)/7r^ with 7r;(xcq) = Pcq/Ph eveutually yielding 

..^HMi^. (9) 

Pip 

Note that the parameters rio and n^ depend on ps, which is consistent with the results 
presented in the right panel of Fig. [31 Indeed, we find both in the numerical simulations 
and in the analytical model that ps has an influence on rent fluctuations, but not on the 
mean rent (which in the analytical model is independent of Uq and n^,). 

To determine the distribution of x, we first write the master equation governing 
the probability Pn{t) that the random walk is at a position n = x/a at time t, using the 
probability rates q and q{x) given in Eq. ([2]): 

dPn{t) 



dt 



-{q + qn)Pn{t) + qPn-noUrit) + gn+n.,n,-Pn+n.,n, (i), (10) 
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with the notation g„ = q{na). If a is small, the distribution essentially appears as 
continuous, and Pn{t) can be described by a continuous distribution (p{x,t) through 
the relation Pn(t) = aip{na,t). Expanding the different terms to second order in a, we 
obtain the Fokker-Planck equation 
d(p{x,t) dJ{x,t) 



dt dx 

with 



11^ 



J{x, t) = alnoUrq — n^niq{x)]ip{x, t) 

~Ydx V-'^^°'^^^^ + nln^Q{x))^{x, t)J . (12) 

Searching for a stationary solution with zero flux, we solve the equation J{x, t) = [25] 
which yields, to leading order in a, a stationary distribution of the form 

(^(x) oce^3(^), (13) 

where g{x) is a function that can be determined explicitly, and that exhibits a maximum 
for Xeq- Given that a is small, (p{x) is dominated by the vicinity of the maximum of 
g{x) and one can expand g{x) around x^q, to second order in x — x^q, leading to 



^{x) = — T^exp 
aVZn 



\X XeqJ 



2a2 



(14) 



with 



(T=i/— ; (noTlr + 12^111). (15) 

V 2 In 10^ ^ ^ ' 

Coming back to the rent variable p, we thus recover a lognormal law, similarly to the 
results obtained from numerical simulations. Fig. [5] shows that this simple analytical 
model is in good agreement with the numerical simulations as far as the equilibrium rent 
is concerned, which is a non-trivial result. The analytically predicted standard deviation 
presents a qualitatively correct behavior as a function of density (namely, it decreases 
when density is increased), but its magnitude is overestimated by a factor of the order of 
3. This discrepancy may be partly explained by the fact that the random walk performed 
by the logarithm of the rent in the numerical simulations is anticorrelated, meaning that 
a jump of characteristic length UoUr to the right is followed by a jump of characteristic 
length nyUi to the left, and reciprocally. This is not the case with the random walk 
we consider here, in which the effective jumps of length Uq and n„ are uncorrelated, 
allowing several consecutive jumps in the same direction. This decorrelation is expected 
to make the rent dispersion larger than in the simulated model. 

An interesting test of this model would be to fit empirical data of rent distributions 
looking like lognormal distributions. Focusing for simplicity on our analytical model, 
measurement oixeq and a on such empirical data would determine respectively TTrPirir/rii 
and apsUi/pi, supposing that the vacancy rate fixes the density p. Some additional 
knowledge of the rental housing market and its dynamics is then needed to find 
acceptable values of the characteristic rents pi and ps, and of the parameters describing 
the "microscopic" dynamics of rent: n^, ni, vr^ and a. Such microscopic parameters may 
however be difficult to determine. 
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5. Conclusion and perspectives 

We have studied a simple model of urban housing market, focusing on the distribution of 
rents (or rents). The mechanisms of rent formation are simple and intuitive: flats which 
are highly demanded have an increasing rent, whereas flats for which agents have a low 
interest have a decreasing rent. The balance between supply and demand resulting in the 
equilibrium rent distribution can be seen as a bargaining between tenants and landlords, 
where very few hypotheses are introduced concerning the respective bargaining powers 
of both categories of agents. From a physics viewpoint, one can see that non-linearities 
encoded in behavioral rules (through the fact that the probability to change flat for 
tenants, or to lower the rent for landlords, depends on the rent) play a major role in 
allowing for a saturation of rent dispersion. Without these non-linearities, the rent 
distribution would continuously spread, just like in the case of an unbounded random 
walk. 

The main specificity of this work when compared to economic literature is the simple 
stochastic behavior of agents, whereas most models in economics use rules which are at 
least partly deterministic. As a result of the use of stochastic behavior rules for agents, 
the equilibrium of our model is different from the equilibria commonly encountered 
in the economic litterature, which are usually static in some sense: once rents have 
reached an equilibrium distribution, firms do not have an incentive to change them, so 
that they do not evolve anymore [TU [8]. In our model, the equilibrium distribution of 
rents is stationary even though rents still evolve at equilibrium. However, the system 
has reached a stationary state, where moves of tenants and rent adjusments by landlords 
are such that they do not change the rent distribution, a result which is reproduced by 
our analytical approach. 

This work also emphasizes the importance of the density of agents, which plays a key 
role in the determination of both the equilibrium rent and the amplitude of fluctuations. 
The equilibrium rent increases with density, while the relative amplitude of fluctuations 
(given by the standard deviation of the logarithm of the rent) decreases with density. 
This tends to validate an hypothesis made in an earlier work [26], of a direct positive 
relationship between equilibrium rent and density. Numerical results are qualitatively 
well reproduced by a model of biased random walk with uncorrelated steps, providing a 
simple understanding of the dynamics of the simulated model. A quantitative agreement 
is even found for the equilibrium rent. 

Extensions of this work could be considered in several directions. Our model could 
be considered as a labor market model with on-the-job search (see [12] )• There is however 
no equivalent of frictional unemployment: the vacancy rate is non zero, but all agents 
have a home, which corresponds to all workers being employed on the labor market. A 
further research direction could be to study the labor market by modifying our model 
to allow for unemployment (that is, agents leaving their jobs even before having found 
another one). Such a work could attempt to explain why wage dispersion is larger in 
real data than predicted by search models [13]. More immediately, introducing new 
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ingredients in the model would make it slightly more realistic. For instance, one could 
introduce a cost of move Ap: agents would move into a new flat only ii p^—p > Ap. This 
cost of move should have a friction-like effect on the housing market, and is expected to 
increase the width of the distribution of rents. In addition, differentiation between flats 
could be introduced, which might allow for a more precise comparison with the literature 
on urban economics, where fiats are differentiated for instance by their distance to the 
city center [20] . 

On the statistical physics side, it would be interesting to investigate in more 
details the connection with ratchet models and dichotomous Markov noise [19]. Indeed, 
the dynamics of the rent of a given fiat switches between exponential increase and 
decrease phases, depending on whether the flat is occupied or empty. This situation 
is similar in spirit to ratchet models with dichotomous Markov noise, in which the 
dynamics of a particle randomly switches between two different regimes, corresponding 
for instance to two different potential energy profiles [19]. There are however some 
subtleties in our model which make it different from a standard ratchet model. First 
of all, the change of state (occupied or empty) of a given flat is not strictly speaking a 
random noise, but results from the A^-body dynamics of the whole system. However, 
this first difficulty could be circumvented by considering a self-consistent noise in a 
mean- field spirit. More importantly, the probability that an agent changes flat depends 
on the corresponding rent, which introduces a non-trivial coupling between rents and 
occupancy. Hence standard methods from ratchet theory cannot be straightforwardly 
applied, and extensions would be necessary. 
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